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Abstract 

An angular analysis and a measurement of the differential branching fraction of 
the decay B® —> 4>p + are presented, using data corresponding to an integrated 
luminosity of 3.0fb _1 of pp collisions recorded by the LHCb experiment at \fs = 7 
and 8TeV. Measurements are reported as a function of q 2 , the square of the dimuon 
invariant mass and results of the angular analysis are found to be consistent with 
the Standard Model. In the range 1 < q 2 < 6GeV 2 /c 4 , where precise theoretical 
calculations are available, the differential branching fraction is found to be more 
than 3 a below the Standard Model predictions. 
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1 Introduction 


The decay B° s —» <pp + p~ is mediated by a b —> s flavour changing neutral current (FCNC) 
transition. In the Standard Model (SM) it is forbidden at tree-level and proceeds via loop 
diagrams as shown in Fig. [l] In extensions of the SM, new heavy particles can appear in 
competing diagrams and affect both the branching fraction of the decay and the angular 
distributions of the final-state particles. 

This decay channel was first observed and studied by the CDF collaboration 01 
and subsequently studied by the LFICb collaboration using data collected during 2011, 
corresponding to an integrated luminosity of 1.0 fb” 1 |3|. While the angular distributions 
were found to be in good agreement with SM expectations, the measured branching fraction 
differs from the recently updated SM prediction by 3.1 a |4|[5]. A similar trend is also seen 
for the branching fractions of other b —> sp + p~ processes, which tend to be lower than 
SM predictions (6]-[8|. 

This paper presents an updated analysis of the decay B J? —> </>(—)■ K + K~)p + pT using 
data accumulated by LHCb in pp collisions, corresponding to an integrated luminosity of 
1.0 fb” 1 collected during 2011 at 7TeV and 2.01b” 1 collected during 2012 at 8TeV centre- 
of-mass energy. The differential branching fraction dB(B® —s- 0/x + /i”)/d q 2 is determined as 
a function of g 2 , the square of the dimuon invariant mass. In addition, a three-dimensional 
angular analysis in cos#;, cos 9k and $ is performed in bins of g 2 . Here, the angle 9k (#z) 
denotes the angle of the K~ (/i”) with respect to the direction of flight of the B° s meson 
in the A' + A'” (p + p~) centre-of-mass frame, and <F denotes the angle between the /i + /i” 
and the K + K~ decay planes in the B° s meson centre-of-mass frame. Compared to the 
previously published fit of the one-dimensional projections of the decay angles |3|, the full 
three-dimensional angular fit gives improved sensitivity and allows access to more angular 
observables. 

The decay B° s —>■ <pp + p~ is closely related to the decay B° —> K*°p + p~, which has 
been studied extensively by LHCb |6,9,10 . Although B ° meson production is suppressed 
with respect to the B° meson by the fragmentation fraction ratio f s /fd ~ 1/4, the 
narrow <p resonance allows a clean selection with low background levels. Furthermore, 
the contribution from the S wave, where the K + K~ system is in a spin-0 configuration, 
is expected to be low 0 Since the K + K p + p final state is not flavour-specific, the 
angular observables accessible in the decay B ° —> <pp + p~ are the CP averages A L , 
and the CP asymmetries ^ 6 , 8,9 12 . The flavour-averaged differential decay rate, as a 
function of the decay angles in bins of q 2 , is given by 


1 d 3 T 9 

dr /dg 2 dcos#; dcos#A' d<F 327T 


[f (1 - F l ) sin 2 6 k + F h cos 2 6 k 


+ | (1 — Fl) sin 2 0 K cos 2 6\ — F L cos 2 d K cos 2 6\ 

+ S 3 sin 2 6 k sin 2 #; cos 2<F + S 4 sin 29 K sin 2#; cos <F 
+ A 5 sin 29 K sin 9i cos <F + A 6 sin 2 9 K cos #; 

+ SV sin 29k sin #; sin <f> + Ag sin 29k sin 2#; sin $ 


( 1 ) 
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Figure 1: Examples of b —> s loop diagrams contributing to the decay B® —> (f)p, + p, in the SM. 


The T-odd CP asymmetries A 8 and Ag are predicted to be close to zero in the SM and 
are of particular interest, as they can be large in the presence of contributions beyond the 
SM [E 


2 Detector and simulation 


The LHCb detector 113141 is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < rj < 5, designed for the study of particles containing b or c quarks. The 
detector includes a high-precision tracking system consisting of a silicon-strip vertex 
detector surrounding the pp interaction region, a large-area silicon-strip detector located 
upstream of a dipole magnet with a bending power of about 4 Tin, and three stations 
of silicon-strip detectors and straw drift tubes placed downstream of the magnet. The 
tracking system provides a measurement of momentum, p, of charged particles with a 
relative uncertainty that varies from 0.5% at low momentum to 1.0% at 200GeV/c. The 
minimum distance of a track to a primary vertex, the impact parameter (IP), is measured 
with a resolution of (15 + 29/px) pm, where px is the component of the momentum 
transverse to the beam, in GeV/c. Different types of charged hadrons are distinguished 
using information from two ring-imaging Cherenkov detectors. Photons, electrons and 
hadrons are identified by a calorimeter system consisting of scintillating-pad and preshower 
detectors, an electromagnetic calorimeter and a hadronic calorimeter. Muons are identified 
by a system composed of alternating layers of iron and multiwire proportional chambers. 
The online event selection is performed by a trigger 1151, which consists of a hardware 
stage, based on information from the calorimeter and muon systems, followed by a software 
stage, which applies a full event reconstruction. 

Simulated signal samples are used to determine the effect of the detector geometry, 
trigger, reconstruction and selection on the signal efficiency. In addition, simulated 
background samples are used to determine the pollution from specific background processes. 


In the simulation, pp collisions are generated using Pythia 16,17 with a specific LHCb 
configuration 18 . Decays of hadronic particles are described by EvtGen 19 , in which 


final-state radiation is generated using PHOTOS 120 1 . The interaction of the generated 
particles with the detector, and its response, are implemented using the Geant 4 toolkit 21 
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22 as described in Ref. 23 . Data-driven corrections are applied to the simulated samples 


to account for imperfect modelling of particle identification performance, the B^ meson 
transverse momentum spectrum and B® vertexing quality, as well as track multiplicity. 


3 Selection of signal candidates 

The B Q S —y 0/r + /i~ signal candidates are required to satisfy the hardware trigger requirement, 
which selects muons with p t > 1.48GeV/c in the 7TeV data and p^ > 1.76GeV/c in the 
8 TeV data. In the subsequent software trigger, at least one of the final-state particles is 
required to have both p t > 0.8GeV/c and IP larger than 100 pm with respect to all of 
the primary pp interaction vertices (PVs) in the event. The tracks of two or more of the 
final-state particles are also required to form a vertex that is significantly displaced from 
any PV. 

Signal candidates are accepted if their reconstructed invariant mass is in the range 
5267 < m(K + K~fi + p~) < 5800MeV/c 2 and the invariant mass of the K + K~ system is 
within 12MeV/c 2 of the known 0 mass [24]. The mass resolutions are 19MeV/c 2 for the 
invariant K + K~p + pr mass and 4MeV/c 2 for the K + K~ invariant mass. The final-state 
particles are required to have significant y 2 P with respect to any PV in the event, where y 2 P 
denotes the change in the y 2 of the PV when reconstructed with or without the considered 
track. The four final-state tracks are then fitted to a common vertex which is required 
to be of good quality and significantly displaced from any PV in the event. The signal 
candidate is required to have small y 2 p with respect to a PV in the event. Furthermore, 
the angle 0dira between the reconstructed B° s momentum and the vector connecting the 
PV with the decay vertex is required to be small. 

To further reduce the combinatorial background, a boosted decision tree (BDT) [25] 
using the AdaBoost algorithm 1261 is employed. The BDT is trained with a sample of 
B Q S — > J/0 (—> p + p~)(j)(—> K + K~ ) decays as a signal proxy and events from t-hc upper 
mass sideband, 5567 < m(K + K~p + ) < 5800MeV/c 2 , as a proxy for t-hc background. 
The discriminating variables of the BDT are the y 2 P of the B° s signal candidate and all 
fmal-state tracks, the B° s transverse momentum, the y 2 of the vertex fit (yvtx)> the flight 
distance significance of the signal candidate, and particle identification information for the 
final-state particles. The BDT selection has an efficiency of 96% for signal events with 
a background rejection of 95%. The total efficiency for signal events, including detector 
geometry, trigger and reconstruction effects is 1.1%. 

The reconstructed B Q S mass versus q 2 for signal candidates after the full selection is 
given in Fig. [ 2 ] The signal decay B® — > (pp + p~ is clearly visible as a vertical band. In the 
q 2 region around the J/\b and ^(2S) masses, the tree-level charmonium decays B® —>■ ,//?/; (f> 
and B Q S —)■ -0(2S')0 dominate. The decay B Q S —)■ is used as a control mode throughout 

the analysis. 


3 





m(K + K fl + Hl [GeV/c 2 ] 


Figure 2: Two-dimensional distribution of q 2 versus the invariant mass of the K + K~ /x + qi~ 
system. The signal decay B® —>• <f>p , + is clearly visible within the ±50 MeV/c 2 interval around 
the Bg mass, indicated by the dashed vertical lines. The horizontal lines denote the charmonium 
regions, where the tree-level decays B®—> and Bg—> S)(j) dominate. 


3.1 Backgrounds 


The decays B® — y ,//?/; (f> and B® — y ip(2S)cj) : primarily originating from b —Y ccs tree-level 
processes, are vetoed by rejecting candidates in the q 2 regions 8.0 < q 2 < 11.0GeV 2 /c 4 
and 12.5 < q 2 < 15.0GeV 2 /c 4 . The B° s —y decay can also constitute a peaking 

background if one of the final-state muons is misidentified as a kaon and vice-versa. This 
background is vetoed by rejecting candidates for which the invariant mass of the K ± /j, t 
system, with the kaon reconstructed under the muon mass hypothesis, is within 45 MeV/c 2 


of the known J/^ meson mass 24 , unless the final-state particles fulfil stringent particle 
identification requirements. After the veto is applied, this background contribution is 
found to be negligible. 

The rare baryonic decay A® —y T(1520)(—> pK~)n + /i~ can mimic the signal decay 
if the proton in the final state is misidentified as a kaon. This potential background is 
vetoed by rejecting events with invariant mass close to the known A® baryon mass 24 


where one kaon has the proton mass hypothesis assigned, unless the kaon passes stringent 
particle identification requirements. Assuming a q 2 dependence following Ref. [27 and 
using B(A° b Ap+fi~) = (0.96 ± 0.29) x 10" 6 [28] as an estimate for the unknown 
)• A( 1520)/r + /i~ branching fraction, a yield of 2.0±0.8 A ° b —>• A(I520)fi + p~ background 
events is expected in the signal region, within 50 MeV/c 2 of the known B° s 


mass 


24 


after the veto. The rare decay B° —» can be a peaking background if the 

pion in the final state is reconstructed as a kaon. After suppressing this background 
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Figure 3: Invariant mass distribution for (left) B®—> 4>ij, + /i~ signal decays, integrated over the 
q 2 bins used, and for (right) the control mode B®—> J/ipp in the K + K~ /r _ final state. The 
signal component is given by the solid blue area, the background component by the shaded red 
area. 


using particle identification information, a yield of 1.7 ± 0.4 events is expected in the 
signal region. The background pollution from A® —> 1520)/i + /r~ and B° —> 

decays is neglected in the fit and treated as a systematic uncertainty. Backgrounds from 
semileptonic b —» c(—* s/i"P M )/i + ^ cascade decays and fully hadronic decays such as 
B°, -(■ D~(-+ K+K- 'Ti )7T + , where hadrons are misidentified as muons, are estimated to 
be small and can therefore be neglected. 


4 Differential branching fraction 


Figure [3] shows the K + K~qi + /j,~ invariant mass distribution for the B° s —> 0/r + /i~ signal 
decay integrated over q 2 , as well as for the control mode B ° —> J/ipcj). To determine the 
B®—} signal yields in bins of q 2 , extended maximum likelihood fits are performed. 

The combinatorial background is described by an exponential function, whilst the signal 
component is modelled with the sum of two Gaussian functions with a common mean 
and a radiative power-law tail toward smaller invariant mass values. The parameters 
describing the signal mass shape are determined from a fit to the B® —> J/^tp control mode. 
The q 2 dependence of the signal mass resolution is accounted for by using scale factors, 
which are determined from simulation. The K + K~/jA fi~ invariant mass distributions 
for the signal decay in bins of q 2 are given in Appendix [Aj the yields and corresponding 
uncertainties are listed in Table [lj Integrating over the q 2 bins, the signal yield is found 
to be 432 ± 24. A fit to the control mode B® > J/ipp, which is used for normalisation, 
gives = 62 033 ± 260 decays. 

The differential branching fraction for a given q 2 bin [gT n , g^ ax ] is calculated according 
to 


d B(B° S —)■ 0/x + /i ) 
d q 2 


q 2 . 

imm 


q2 

^max 


Nj/tluf, 


— ■ B(B°^ iV), 


( 2 ) 


where and Nj denote the yield of the signal and normalisation mode, and e^ 
and ej/^^ their respective efficiencies. The branching fractions are given by t3(J/ip —> 
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fi + fi~) = (5.961 ± 0.033) x 10~ 2 [24] and B(B° S ^ J/^4>) = (10.76 ± 0.81) x 10~ 4 . For the 
branching fraction of the normalisation channel B® —>■ J/0 </> the LHCb measurement 


is recalculated using an updated measurement of f s /fd = 0.259 ± 0.015 29 . A weighted 
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average is calculated by combining this updated measurement with the measurements 
by Belle 


30 and CDF 


The resulting relative and absolute differential branching 

fractions are given in Table [lj 

The differential branching fraction is also shown in Fig. [4j overlaid with SM predictions 
from Refs. |4.j5]. In the q 2 region 1.0 < q 2 < 6.0 GeV 2 /c 4 the measured differential branching 
fraction lies 3.3 a below the SM expectation of (4.81 ±0.56) x 10“ 8 GeV“ 2 c' 


32 


For the 

SM predictions, the form factors are determined in a combined fit to the results of light-cone 
sum rule calculations at low q 2 [h| and lattice QCD calculations at high q 2 33 34 . Standard 
Model predictions for the branching fraction at high q 2 that exclusively use the results 
from lattice calculations 35 are found to be larger than the results from the combined 


fit. Owing to their proximity to the charmonium resonances, no predictions are available 
corresponding to the q 2 bins 5.0 < q 2 < 8.0GeV 2 /c 4 and 11.0 < q 2 < 12.5 GeV 2 /c 4 . 

The total branching fraction of the signal decay is given by the integral over the six q 2 
bins. To account for the fraction of signal events in the vetoed q 2 regions, a correction 
factor / veto = 1.520 ± 0.003 ± 0.043 is applied, which is determined using the calculation in 
Ref. [36| with updated form factors from Ref. [37|. The first given uncertainty is statistical, 
the second is systematic. 

The resulting relative and total branching fractions are 


B(B^ 0/i + /i ) 
B(B° S -> M4>) 
B(B° S -4 <MV) 


(7.41+°; 42 ± 0.20 ± 0.21) x 10~ 4 , 
(7.971^3 ± 0.22 ± 0.23 ± 0.60) x 10“ 7 , 


where the uncertainties are (from left to right) statistical, systematic, and from the 
extrapolation to the full q 2 region. For the total branching fraction, a further uncertainty 
originates from the uncertainty on the branching fraction of the normalisation mode. 


4.1 Systematic uncertainties 

For the branching fraction ratio £>(£>]?—)• (ftn + p~)/B(B® —* J/V>0), systematic uncertainties 
are mostly due to uncertainties on the efficiency ratio which is taken from 

simulation. To evaluate the size of uncertainties affecting the efficiency ratio, it is recal¬ 
culated after applying the corresponding systematic variation to the simulated samples. 
The observed deviation is taken as systematic uncertainty. The procedure to correct the 
tracking efficiency in simulation introduces a systematic uncertainty on the efficiency ratio 
of less than 0.6%. The correction to particle identification performance in simulation 
has a systematic uncertainty of 0.5%. The relative efficiency is further affected by the 
data-driven corrections to the simulation in the distribution of the variables Pt(B®) and 
Xy tx (I?°), as well as the track multiplicity, which have a combined systematic effect of 1.0%. 
The non-uniform angular acceptance detailed in Sec. [5] introduces a dependence of the 
signal efficiency on the underlying physics model. Its effect on the branching fraction 
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Table 1: The signal yields for R° —> 0/r+/i~ decays, as well as the differential branching fraction 
relative to the normalisation mode and the absolute differential branching fraction, in bins of q 2 . 
The given uncertainties are (from left to right) statistical, systematic, and the uncertainty on the 
branching fraction of the normalisation mode. 


q 2 bin [GeV 2 /c 4 ] [10-° GeV~V] 


} [10“ 8 GeV“ 2 c 4 ] 


0.1 < q 2 < 2.0 
2.0 < q 2 < 5.0 60lg° 

5.0 < q 2 < 8.0 83l^ 

11.0 < q 2 < 12.5 70l^ 

15.0 < q 2 < 17.0 83+^ 

17.0 < q 2 < 19.0 54+® 


5.44+°;®® ±0.13 
2.38ig;p ± 0.06 
2.98+g;3g ± 0.07 
4.37l^±0.14 
4.20l^±0.11 
3.68+°;5 q ± 0.13 


5.85+°;gg ± 0.14 ±0.44 
2.56±g;|| ±0.06 ±0.19 
3.2lt°;^± 0.08 ±0.24 
4.71±g;|§ ± 0.15 ±0.36 
4.52+°;®J± 0.12 ±0.34 
3.96^;^ ±0.14 ±0.30 
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Figure 4: Differential branching fraction of the decay F>° —> </>/r+^ _ , overlaid with SM predic¬ 
tions |4|5] indicated by blue shaded boxes. The vetoes excluding the charmonium resonances are 
indicated by grey areas. 


measurement is evaluated by varying the Wilson coefficient Cg used in the generation 
of simulated signal events. By allowing a New Physics contribution of —1.5, which is 
motivated by the global fit results in Ref. 38 , the resulting systematic uncertainty is 
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found to be less than 1.6%. The selection requirements introduce a decay-time dependence 
of the efficiencies which can, due to the sizeable lifetime difference in the B® system 39 , 
affect the measured branching fraction [401 The systematic uncertainty is determined with 
simulated B ° —>■ 0/r + /i” signal events, generated using time-dependent decay amplitudes 


as described in Ref. 1121. When varying the Wilson coefficients, the size of the effect is 
found to be at most 1.6%, which is taken as the systematic uncertainty. The statistical 
uncertainty due to the limited size of the simulated signal samples leads to a systematic 
uncertainty of 1.9%. 

The systematic uncertainties due to the parametrisation of the mass shapes are 
evaluated using pseudoexperiments. For the signal mass model, events are generated 
using a double Gaussian mass shape, and then fitted using both the double Gaussian as 
well as the nominal signal mass shape, taking the observed deviation as the systematic 
uncertainty. For the parametrisation of the combinatorial background, the nominal 
exponential function is compared with a linear mass model. The systematic uncertainties 
due to the modelling of the signal and background mass shape are 2.1% and 1.6%, 
respectively. Peaking backgrounds are neglected in the fit for determination of the signal 
yields. The main sources of systematic uncertainty are caused by contributions from 
the decays A® —> pK~fi + /i~ and B° —» K*°p + /d~, resulting in systematic uncertainties of 
0.2 — 2.2%, depending on the q 2 bin. Finally, the uncertainty on the branching fraction of 
the decay J/0 —» /r + /i~ amounts to a systematic uncertainty of 0.6%. The complete list of 
systematic uncertainties is given in Table [2} 

For the total branching fraction of the signal decay, the uncertainty on the branching 
fraction of the normalisation channel is the dominant systematic uncertainty, at the level 
of 7.5%. The uncertainty on the correction factor / veto to account for signal events that 
are rejected by the charmonium vetoes is estimated by varying the Wilson coefficients and 
form-factor parameters, leading to a systematic uncertainty of 2.9%. 


5 Angular analysis 

For the determination of the four CP averages iffi, £ 3 , 4,7 and the four CP asymmetries 
R 5 , 6 , 8,9 an unbinned maximum likelihood fit to the three-dimensional angular distribution 
and the K + K~p + /i _ invariant mass distribution is performed in each q 2 bin. The models 
described in Sec. [4] are used to parametrise the mass line shapes for signal and background. 
The angular distribution of the signal component is given by Eq. [lj The angular background 
distribution is described by the product of second-order Chebyshev polynomials in the 
three decay angles. 

The non-uniform efficiency due to the reconstruction, triggering and selection of signal 
candidates distorts the angular distributions of the final-state particles, as well as the q 2 
distribution. This acceptance effect is parametrised using Legendre polynomials, according 
to 


e(cos 9i, cos 9k, $, q 2 ) = E CklmnPk(cOs9 l )Pi(cOs9 K )P m ($)P n (q 2 ), (3) 

klmn 
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Table 2: Systematic uncertainties [10 5 GeV 2 c 4 ] on the branching fraction ratio d£>(R|? —> 
4>/j, + —> J/i/j(j>)dq 2 per bin of q 2 [GeV 2 /c 4 ]. 


Source 

[0.1,2] 

[2,5] 

[5,8] 

[1142.5] 

[15,17] 

[17,19] 

[1,6] 

[15,19] 

Simulation corr. 

0.01 

0.01 

0.01 

0.01 

0.05 

0.04 

0.00 

0.04 

Angular model 

0.04 

0.00 

0.01 

0.00 

0.01 

0.06 

0.00 

0.01 

Efficiency ratio 

0.06 

0.03 

0.03 

0.06 

0.06 

0.07 

0.02 

0.04 

Signal mass model 

0.02 

0.01 

0.03 

0.03 

0.03 

0.00 

0.05 

0.05 

Bkg. mass model 

0.02 

0.02 

0.02 

0.02 

0.03 

0.05 

0.01 

0.06 

Time acceptance 

0.09 

0.04 

0.05 

0.07 

0.07 

0.06 

0.04 

0.06 

->• p + p“) 

0.03 

0.01 

0.02 

0.02 

0.02 

0.02 

0.01 

0.02 

Peaking bkg. 

0.03 

0.02 

0.02 

0.10 

0.02 

0.01 

0.02 

0.01 

Quadratic sum 

0.13 

0.06 

0.07 

0.14 

0.11 

0.13 

0.07 

0.12 


where Pi(x) denote Legendre polynomials of order i and Cki mn the coefficients that are 
determined by performing a moments analysis using a large sample of simulated B Q S —> 
0/x + /i - signal events generated according to a phase-space model. The maximum order of 
the polynomials that is included is four for cos 9i, two for cos6%, six for the angle $ and 
five for q 2 . In addition, the acceptance is assumed to be symmetric in the decay angles. 
This choice corresponds to the lowest orders of polynomials that describe the acceptance 
effect. The acceptance description is cross-checked using the control mode B° s —>■ J/00. 
An angular analysis of the control mode is performed and the angular observables are 
found to be in good agreement with the previous measurement 
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Appendix [B] gives the one-dimensional angular distributions of the signal decay in 
each q 2 bin, overlaid with the projections of the likelihood fit. For the q 2 bins with the 
lowest number of signal candidates, pseudoexperiments show the likelihood estimator to 
be biased for certain observables due to boundary effects which arise from the requirement 
of Eq. [l] being positive for all values of the decay angles. Therefore, the Feldman-Cousins 
method [411 is used to determine confidence regions for each observable, which guarantees 
correct coverage for low signal yields. The remaining signal observables are treated as 
nuisance parameters following the plugin method [42|. The Feldman-Cousins scans for the 


angular observables in bins of q 2 are given in Appendix [cj In some q 2 bins, the shape of 
the obtained confidence level is dominated by effects arising from the boundary condition. 
Table [3] gives the minima of the Feldman-Cousins scans and the 68% confidence intervals. 
The linear correlations between the angular observables in the different q 2 bins are given 
D The angular observables are shown in Fig. [5j overlaid with SM predictions 
5 . No predictions are given for S 7 and A 5)6i8!9 ; they are expected to be close 


in Appendix 
from Refs. 14 


to zero in the SM. 
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Table 3: (Top) CP-averaged angular observables Fl and 63 , 4,7 and (bottom) CP asymmetries 
^ 5 , 6 , 8,9 obtained from the unbinned maximum likelihood fit, where the first uncertainty is 
statistical and the second is systematic. 


q 2 bin 

[ GeV 2 / c 4 ] 

Fl 



63 



64 



67 



0.1 

< 

q 2 

< 

2.0 

0-20 18:88 

± 

0.02 

- 0.051818 

± 

0.01 

n 97+O.28 
u - z '- 0.18 

± 

0.01 

0 . 0418:18 

± 

0.00 

2.0 

< 

q 2 

< 

5.0 

0-68±g;il 

± 

0.03 

-o-o6i8i8 

± 

0.01 

—0 47+ 030 

u -^'- 0.44 

± 

0.01 

- 0.031818 

± 

0.01 

5.0 

< 

q 2 

< 

8.0 

0 . 5418 $ 

± 

0.02 

—0 1 n+°- 20 

u - iu - 0.29 

± 

0.01 

—0 1 o+°- 15 
U - 1 U - 0.18 

± 

0.01 

0.041818 

± 

0.01 

11.0 

< 

q 2 

< 

12.5 

0 - 2918:1 

± 

0.04 

— 0 1Q+ 0 - 20 
u - la - 0.23 

± 

0.01 

—0 47+ 0 21 

u -^'- 0.29 

± 

0.01 

0.00I8I® 

± 

0.01 

15.0 

< 

q 2 

< 

17.0 

f) oq+ 0-09 
u - zo - 0.08 

± 

0.02 

— 0 flf>+ 0 - 16 
u - uo - 0.19 

± 

0.01 

—0 04+°- 15 
u - uo - 0.15 

± 

0.01 

0 19+°- 16 
u - iz - 0.13 

± 

0.01 

17.0 

< 

q 2 

< 

19.0 

0 40 +0 ' 13 
U - 4 U - 0.15 

± 

0.02 

—0 07+°- 23 
U.U(_o 27 

± 

0.02 

— 0 3 Q+ 0 - 25 
U - OJ - 0.34 

± 

0.02 

0 ?0 +a29 
u - zu -0.22 

± 

0.01 

1.0 

< 

q 2 

< 

6.0 

0 fi 3 + 0 - 09 
u - oo - 0.09 

± 

0.03 

-o.o 218;18 

± 

0.01 

— 0 1 Q+ 0 - 14 
u - la - 0.13 

± 

0.01 

—0 04+°- 14 
u - uo - 0.14 

± 

0.00 

15.0 

< 

q 2 

< 

19.0 

Q OQ+ 0 - 07 
u - za - 0.06 

± 

0.02 

- 0 . 0918:11 

± 

0.01 

—fl 14+0+1 
u - 14 -0.11 

± 

0.01 

fl 1 Q+0.11 

u - 10 -o.ii 

± 

0.01 


q 2 bin 

[GeV : 

! /c 4 ] 

A 



Ae 



4^8 
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0.1 

< 

q 2 

< 

2.0 

—0 02 +013 
u - uz - 0.13 

± 

0.00 

- 0 . 1918:18 

± 

0.01 

n i n+°- 14 

u - iu - 0.14 

± 

0.00 

0 . 0318:11 

± 

0.01 

2.0 

< 

q 2 

< 

5.0 

0.091818 

± 

0.01 

0.091818 

± 

0.02 

fl iq+ 0.26 
u - ia -0.21 

± 

0.01 

- 0-131818 

± 

0.01 

5.0 

< 

q 2 

< 

8.0 

0.041811 

± 

0.01 

—fl ni+ 0 - 14 
U ' U1 -0.12 

± 

0.01 

—0 12 +al7 
u - iz - 0.19 

± 

0.01 

—0 04+°- 17 
u - uo - 0.16 

± 

0.01 

11.0 

< 

q 2 

< 

12.5 

0.081811 

± 

0.01 

— fl 1fi+°- 16 
u - lo - 0.18 

± 

0.01 

-0.01I8I® 

± 

0.01 

—0 02 + °' 16 
u - uz - 0.15 

± 

0.01 

15.0 

< 

q 2 

< 

17.0 

o.o2i8:l 3 

± 

0.01 

fl fll +°- 12 
U - Ui - 0.17 

± 

0.01 

0.0818:1® 

± 

0.01 

fl 91 + 0.18 
U - Zi -0.12 

± 

0.01 

17.0 

< 

q 2 

< 

19.0 

n 1 q+ 0.29 
u ' iO - 0.27 

± 

0.01 

- 0 . 0418:18 

± 

0.01 

-0.16I8I8 

± 

0.01 

—fl 02 + °' 19 
u - uz - 0.19 

± 

0.01 

1.0 

< 

q 2 

< 

6.0 

0 2fl +0 ' 13 
u - zu - 0.13 

± 

0.00 

0.08181? 

± 

0.01 

-0.00I8I® 

± 

0.00 

—0 fli+°- 13 
U - U 1 - 0.13 

± 

0.01 

15.0 

< 

q 2 

< 

19.0 

fl 1 1 +0.10 
u - 4i -0.10 

± 

0.00 

0.00I8I? 

± 

0.01 

0.031818 

± 

0.00 

fl 1 o+O.ll 
u - lz - 0.09 

± 

0.00 


5.1 Systematic uncertainties 

The systematic uncertainties on the angular observables are evaluated using large numbers 
of pseudoexperiments where simulated events are generated to reflect the measured angular 
distributions and event yield. To reduce statistical effects, each sample is fitted twice, once 
with and once without systematic variations. 

The angular acceptance correction, which is determined from simulation, is a significant 
source of systematic uncertainties for the angular observables. The data-driven correc¬ 
tions of the distributions of yy tx (.B?) an d track multiplicity, as well as particle 

identification performance and tracking efficiency, amount to a systematic uncertainty 
of less than 0.01 in total. Furthermore, the kinematic distributions of the final-state 
particles are cross-checked using the control mode B Q S —y J/xf: (f). Correcting for kinematic 
differences between data and simulation amounts to a systematic deviation of less than 
0.01. The effect of the limited size of the simulated signal samples on the acceptance 
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Figure 5: CP-averaged angular observables Pl and 53,4,7 and CP asymmetries 6,8 ,9 shown by 
black dots, overlaid with SM predictions [4,5], where available, indicated as blue shaded boxes. 
The vetoes excluding the charmonium resonances are indicated by grey areas. 
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description is evaluated by varying the Legendre coefficients in Eq. [3] according to their 
corresponding covariance matrix. The resulting systematic uncertainty is smaller than 0.02 
for all observables and q 2 bins. The four-dimensional acceptance correction is evaluated 
at the centre of each q 2 bin. To estimate the systematic effect due to this, an alternative 
acceptance description is used, where a separate three-dimensional acceptance is used for 
each q 2 bin. The resulting systematic deviation is negligible. 

The systematic effect of neglecting peaking backgrounds is evaluated by performing 
toy studies, where simulated —> pK~ p + pT and B° — » K*°p + p~ background events are 
added, according to their expected yields for the specific q 2 bin. The resulting systematic 
deviations of the angular observables are smaller than 0.01 for all observables and q 2 bins. 


The S-wave pollution for the decay B° s —y 0/L 
the B®—> decay, at the level of 1.1% 


11 


/i - is expected to be similar to that of 
The effect of the S wave in the K + K~ 
system on the angular observables is determined to be smaller than 0.01 using toy studies. 
The combinatorial background is described using second-order Chebyshev polynomials 
determined from the upper mass sideband. The systematic uncertainty associated with 
this model choice is estimated by using first-order polynomials as an alternative. With 
a systematic effect of up to 0.04 on the angular observables, depending on q 2 bin, this 
constitutes the dominant systematic uncertainty for the angular analysis. In addition, the 
effect of fixing the angular background parameters in the nominal fit is evaluated using 
toy studies. The systematic deviation is found to be smaller than 0.02 for all observables 
and q 2 bins. 

The total systematic uncertainty, given by the quadratic sum over all systematic effects, 
is found to be small compared to the statistical uncertainties for all angular observables in 
all q 2 bins. 


6 Conclusions 

Measurements of the differential branching fraction and the first full three-dimensional 
angular analysis of the decay B ° —> (j)p + are presented, using data collected by the 
LHCb experiment in pp collisions, corresponding to an integrated luminosity of 3.0fb _1 . 
The results are given in Tables [I] and [3] and are the most precise measurements of these 
quantities to date. The CP-averaged angular observables S 4 and S 7 are determined for the 
first time for this decay. The determination of the CP asymmetries A 5 and A s constitutes 
the first measurement of these quantities for any rare b —>• sil decay, providing additional 
constraints in global fits. All angular observables are found to be compatible with SM 
predictions. 

The Pj? — y 0/P pr branching fraction relative to the normalisation mode B® — y J/ij. ; 0 is 
measured to be 


P(Pg->0/P/i-) 
B{B. J/00) 


(7Alt° 0 il ± 0.20 ± 0 . 21 ) x 10 “ 4 , 
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and the resulting total absolute branching fraction is measured to be 


= (7.97±g;tl ± 0.22 ± 0.23 ± 0.60) x 10 


-7 


where the uncertainties are (from left to right) statistical, systematic, and from the 
extrapolation to the full q 2 region. For the total branching fraction, a further uncertainty 
originates from the uncertainty on the branching fraction of the normalisation mode. 
The measured branching fraction is compatible with the previous measurement |3| and 
lies below SM expectations. For the q 2 region 1.0 < q 2 < 6.0GeV 2 /c 4 the differential 
branching fraction of (2.58lg;3f ± 0.08 ± 0.19) x 10 
SM prediction of (4.81 ± 0.56) x 10 -8 GeV~ 2 c 4 |4]|5 1 32 . 


GeV 2 c 4 is more than 3 a below the 
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Appendices 

A Invariant mass distributions 
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Figure 6: Invariant mass distributions for B® —> 0/i + /r signal decays in bins of q 2 . The signal 
component is shown by the solid blue area, the background component by the shaded red area. 
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B Angular fit projections 
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Figure 7: One-dimensional projections of the fit to the angles cos 6 i, cos 9k, in bins of q 2 . The 
signal component is shown by the solid blue area, the background component by the shaded red 


area. 
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Figure 8: One-dimensional projections of the fit to the angles cos 6 i, cos 9k, in bins of q 2 . The 
signal component is shown by the solid blue area, the background component by the shaded red 


area. 
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C Confidence intervals 



Figure 9: Confidence level obtained from a likelihood scan (shaded blue) and from the Feldman- 
Cousins method (solid black). The shaded blue and solid red vertical lines indicate the corre¬ 
sponding 68% CL intervals obtained from the likelihood scan and the Feldman-Cousins method, 
respectively. 
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Figure 10: Confidence level obtained from a likelihood scan (shaded blue) and from a Feldman- 
Cousins method (solid black). The shaded blue and solid red vertical lines indicate the corre¬ 
sponding 68% CL intervals obtained from the likelihood scan and the Feldman-Cousins method, 
respectively. 
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Figure 11: Confidence level obtained from a likelihood scan (shaded blue) and from a Feldman- 
Cousins method (solid black). The shaded blue and solid red vertical lines indicate the corre¬ 
sponding 68% CL intervals obtained from the likelihood scan and the Feldman-Cousins method, 
respectively. 
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Figure 12: Confidence level obtained from a likelihood scan (shaded blue) and from a Feldman- 
Cousins method (solid black). The shaded blue and solid red vertical lines indicate the corre¬ 
sponding 68% CL intervals obtained from the likelihood scan and the Feldman-Cousins method, 
respectively. 
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Figure 13: Confidence level obtained from a likelihood scan (shaded blue) and from a Feldman- 
Cousins method (solid black). The shaded blue and solid red vertical lines indicate the corre¬ 
sponding 68% CL intervals obtained from the likelihood scan and the Feldman-Cousins method, 
respectively. 
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Figure 14: Confidence level obtained from a likelihood scan (shaded blue) and from a Feldman- 
Cousins method (solid black). The shaded blue and solid red vertical lines indicate the corre¬ 
sponding 68% CL intervals obtained from the likelihood scan and the Feldman-Cousins method, 
respectively. 
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Figure 15: Confidence level obtained from a likelihood scan (shaded blue) and from a Feldman- 
Cousins method (solid black). The shaded blue and solid red vertical lines indicate the corre¬ 
sponding 68% CL intervals obtained from the likelihood scan and the Feldman-Cousins method, 
respectively. 
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Figure 16: Confidence level obtained from a likelihood scan (shaded blue) and from a Feldman- 
Cousins method (solid black). The shaded blue and solid red vertical lines indicate the corre¬ 
sponding 68% CL intervals obtained from the likelihood scan and the Feldman-Cousins method, 
respectively. 
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D Correlation matrices 


Table 4: Correlation matrices for the q 2 bins 0.1 < q 2 < 2.0GeV 2 /c 4 , 2.0 < q 2 < 5.0GeV 2 /c 4 and 
5.0 <q 2 < 8.0 GeV 2 /c 4 . 


Correlation matrix for 0.1 < q 2 < 2.0 GeV 2 /c 4 



F l 

S3 

s 4 

A5 

Aq 

s 7 

As 

Ag 

Fl 

1.00 

0.03 

-0.15 

— 

0.02 

0.10 

0.03 

- 

s 3 


1.00 

0.04 

0.07 

— 

0.05 

-0.18 

-0.05 

S 4 



1.00 

-0.13 

-0.09 

-0.19 

0.06 

-0.09 

A5 




1.00 

0.11 

0.06 

-0.14 

0.10 

Ag 





1.00 

0.07 

-0.03 

-0.16 

s 7 






1.00 

-0.30 

0.03 

As 







1.00 

0.06 

Ag 








1.00 


Correlation matrix for 2.0 < q 2 < 5.0GeV 2 /c 4 



Fl 

S3 

S 4 

A5 

Ag 

Sr 

^8 

Ag 

Fl 

1.00 

-0.05 

0.27 

0.04 

-0.09 

0.02 

0.02 

-0.16 

S3 


1.00 

-0.23 

-0.06 

-0.05 

0.20 

-0.11 

0.40 

s 4 



1.00 

0.11 

0.16 

0.14 

-0.41 

-0.33 

A5 




1.00 

-0.24 

-0.31 

0.06 

0.08 

A& 





1.00 

-0.03 

0.05 

0.11 

s 7 






1.00 

-0.05 

-0.02 

^8 







1.00 

-0.16 

Ag 








1.00 


Correlation matrix for 5.0 < q 2 < 8.0GeV 2 /c 4 



Fl 

S3 

s 4 

A5 

^6 

Sr 

^8 

Ag 

Fl 

1.00 

- 0.03 

- 0.01 

0.11 

0.01 

0.03 

- 0.07 

0.09 

S3 


1.00 

- 0.03 

- 0.01 

- 0.02 

- 0.18 

— 

0.05 

S 4 



1.00 

- 0.05 

- 0.01 

- 0.03 

0.20 

- 0.08 

A5 




1.00 

— 

0.14 

- 0.05 

- 0.16 

Ag 





1.00 

— 

0.05 

- 0.25 

Sr 






1.00 

0.04 

- 0.14 

As 







1.00 

- 0.04 

Ag 








1.00 
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Table 5: Correlation matrices for the q 2 bins 11.0 < q 2 < 12.5GeV 2 /c 4 , 15.0 < q 2 < 17.0GeV 2 /c 4 
and 17.0 < q 2 < 19.0 GeV 2 /c 4 . 


Correlation matrix for 11.0 < q 2 < 12.5GeV 2 /c 4 



Fl 

S3 

5 4 

A5 

^6 

S7 

^8 

A9 

Fl 

1.00 

0.21 

0.11 

0.20 

0.06 

0.21 

-0.05 

0.03 

Ss 


1.00 

0.02 

0.12 

0.18 

0.08 

-0.02 

0.04 

5 4 



1.00 

-0.11 

-0.37 

0.26 

-0.01 

-0.09 

M 




1.00 

-0.16 

-0.07 

0.22 

0.04 

^6 





1.00 

-0.04 

-0.11 

0.14 

S7 






1.00 

-0.36 

0.08 

As 







1.00 

-0.23 

A9 








1.00 


Correlation matrix for 15.0 < q 2 < 17.0 GeV 2 /c 4 



F l 

S 3 

54 

A 5 

^6 

S 7 

^8 

A 9 

Fl 

1.00 

-0.01 

- 0.01 

0.07 

-0.05 

0.06 

-0.06 

-0.06 

S 3 


1.00 

0.03 

-0.06 

0.11 

-0.08 

-0.06 

0.15 

S 4 



1.00 

0.01 

-0.07 

0.04 

0.22 

0.04 

A 5 




1.00 

-0.04 

0.14 

0.05 

0.01 

^6 





1.00 

0.05 

0.01 

-0.09 

S 7 






1.00 

-0.03 

0.06 

As 







1.00 

-0.11 

A 9 








1.00 


Correlation matrix for 17.0 < q 2 < 19.0 GeV 2 /c 4 



Fl 

5 3 

5 4 

A 5 

^6 

5 7 

^8 

A 9 

F l 

1.00 

-0.12 

0.08 

-0.32 

0.06 

-0.04 

0.15 

-0.01 

S 3 


1.00 

-0.04 

0.40 

0.16 

0.22 

— 

0.02 

5 4 



1.00 

0.25 

0.30 

-0.13 

0.40 

0.14 

A 5 




1.00 

-0.05 

0.16 

-0.06 

0.28 

Aq 





1.00 

-0.03 

0.19 

-0.05 

S 7 






1.00 

-0.02 

0.18 

As 







1.00 

-0.02 

A 9 








1.00 
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Table 6: Correlation matrices for the q 2 bins 1.0 < q 2 < 6.0 GeV 2 /c 4 and 15.0 < q 2 < 19.0 GeV 2 /c 4 . 


Correlation matrix for 1.0 < q 2 < 6.0GeV 2 /c 4 



Fl 

S3 

s 4 

A5 

^6 

£7 

As 

^9 

Fl 

1.00 

-0.02 

0.08 

0.08 

-0.04 

0.07 

0.03 

-0.05 

s 3 


1.00 

-0.07 

0.13 

-0.13 

0.10 

-0.06 

0.12 

S 4 



1.00 

0.17 

0.10 

-0.06 

-0.01 

-0.13 

A5 




1.00 

-0.14 

-0.05 

-0.05 

0.04 

Ag 





1.00 

-0.04 

-0.09 

-0.06 

S 7 






1.00 

0.07 

0.04 

^8 







1.00 

-0.10 

Ag 








1.00 


Correlation matrix for 15.0 < q 2 < 19.0 GeV 2 /c 4 



F l 

S3 

s 4 

A5 

Ag 

S 7 

^8 

Ag 

Fl 

1.00 

0.01 

-0.05 

-0.03 

— 

0.07 

-0.04 

-0.08 

S3 


1.00 

-0.03 

0.16 

0.09 

0.05 

0.02 

0.07 

s 4 



1.00 

0.12 

0.03 

0.03 

0.24 

0.09 

A5 




1.00 

-0.09 

0.08 

0.04 

0.14 

Ag 





1.00 

0.05 

0.05 

-0.11 

S7 






1.00 

0.01 

0.10 

As 







1.00 

-0.07 

A9 








1.00 


27 



References 

[1] CDF collaboration, T. Aaltonen et al., Measurement of the forward-backward asym¬ 
metry in the B —» K k)/i + /i~ decay and first observation of the B® —» decay, 

Phys. Rev. Lett. 106 (2011) 161801, arXiv: 1101.1028. 

[2] CDF collaboration, T. Aaltonen et al., Observation of the baryonic flavor-changing 
neutral current decay A & —» A/i + /i”, Phys. Rev. Lett. 107 (2011) 201802, 
arXiv:1107.3753, 

[3] LHCb collaboration, R. Aaij et al. , Differential branching fraction and angular analysis 
of the decay B° s —>■ f>p + , JHEP 07 (2013) 084, arXiv: 1305.2168. 

[4] W. Altmannshofer and D. M. Straub, Aere physics in b —» s transitions after LHC 
run 1 , Eur. Phys. J. C75 (2015) 382, arXiv: 1411.3161. 

[5] A. Bharucha, D. M. Straub, and R. Zwicky, R — y V£ + £~ in the Standard Model from 
light-cone sum rules, arXiv: 1503.05534. 

[6] LHCb collaboration, R. Aaij et al, Differential branching fraction and angular analysis 
of the decay R° —» K*°n + n~, JHEP 08 (2013) 131, arXiv: 1304.6325. 

[7] LHCb collaboration, R. Aaij et al., Differential branching fractions and isospin 

asymmetries of B —>• decays , JHEP 06 (2014) 133, arXiv: 1403.8044. 

[8] LHCb collaboration, R. Aaij et al, First observations of the rare decays B + —* 
K + 7T + 7r~ /a + and B + — > <pK + p + pT , JHEP 10 (2014) 064, arXiv: 1408.1137, 

[9] LHCb collaboration, R. Aaij et al., Measurement of form-factor-independent, ob¬ 
servables in the decay R° —> K* 0 p + p~, Phys. Rev. Lett. Ill (2013) 191801, 
arXiv:1308.1707, 

[10] LHCb collaboration, Angular analysis of the B° —> K*°/j, + /j,~ decay, LHCb-CONF- 
2015-002. 

[11] LHCb collaboration, R. Aaij et al, Amplitude analysis and branching fraction mea¬ 
surement of B° s — > J/fi>K + K~, Phys. Rev. D87 (2013) 072004, arXiv: 1302.1213. 

[12] C. Bobeth, G. Hiller, and G. Piranishvili, CP asymmetries in B —>• /?*(—> Ktt)££ 
and untagged B s , B s —> </>(—)• K + K~)££ decays at NLO, JHEP 07 (2008) 106, 
arXiv:0805.2525. 

[13] LHCb collaboration, A. A. Alves Jr. et al, The LHCb detector at the LHC, JINST 3 
(2008) S08005. 

[14] LHCb collaboration, R. Aaij et al., LHCb detector performance, Int. J. Mod. Phys. 
A30 (2015) 1530022, arXiv: 1412.6352. 


28 








[15] R. Aaij et al, The LHCb trigger and its performance in 2011 , JINST 8 (2013) P04022, 

arXiv:1211.3055. 

[16] T. Sjostrand, S. Mrenna, and P. Skands, PYTHIA 6.f physics and manual, JHEP 05 
(2006) 026, arXiv:hep-ph/0603175, 

[17] T. Sjostrand, S. Mrenna, and P. Skands, A brief introduction to PYTHIA 8.1, Comput, 
Phys. Commun. 178 (2008) 852, arXiv:0710.3820. 

[18] I. Belyaev et al, Handling of the generation of primary events in Gauss, the LHCb 
simulation framework, J. Phys. Conf. Ser. 331 (2011) 032047, 

[19] D. J. Lange, The EvtGen particle decay simulation package, Nucl. Instrum. Meth. 
A462 (2001) 152, 

[20] P. Golonka and Z. Was, PHOTOS Monte Carlo: A precision tool for QED corrections 
in Z and W decays, Eur. Phys. J. C45 (2006) 97, arXiv :hep-ph/0506026. 

[21] Geant4 collaboration, J. Allison et al, Geantf developments and applications, IEEE 
Trans. Nucl. Sci. 53 (2006) 270. 

[22] Geant4 collaboration, S. Agostinelli et al, Geantf: a simulation toolkit, Nucl. Instrum. 
Meth. A506 (2003) 250. 

[23] M. Clemencic et al, The LHCb simidation application, Gauss: design, evolution and 
experience, J. Phys. Conf. Ser. 331 (2011) 032023, 

[24] Particle Data Group, K. A. Olive et al, Review of particle physics, Chin. Phys. C38 
(2014) 090001. 

[25] L. Breiman, J. H. Friedman, R. A. Olshen, and C. J. Stone, Classification and 
regression trees, Wadsworth international group, Belmont, California, USA, 1984. 

[26] R. E. Schapire and Y. Freund, A decision-theoretic generalization of on-line learning 
and an application to boosting, J. Comput. Syst. Sci. 55 (1997) 119. 

[27] L. Mott and W. Roberts, Rare dileptonic decays of A& in a quark model, Int. J. Mod. 
Phys. A27 (2012) 1250016, arXiv: 1108.6129. 

[28] LHCb collaboration, R. Aaij et al, Measurement of the differential branching fraction 
of the decay A° —* Ap + /i~, Phys. Lett. B725 (2013) 25, arXiv: 1306.2577, 

[29] LHCb collaboration, R. Aaij et al, Measurement of the fragmentation fraction 
ratio fs/fd and its dependence on B meson kinematics, JHEP 04 (2013) 001, 
arXiv: 1301.5286, f s /fd value updated in LHCb-CONF-2013-011. 

[30] Belle collaboration, F. Thorne et al, Measurement of the decays B® — > 
J/'0/(lO2O), B° s —* 7 /^/ 2 ( 1525 ) and B° s —> J/ifK + K~ at Belle, Phys. Rev. D88 
(2013) 114006, arXiv:1309.0704. 


29 











[31] CDF collaboration, F. Abe et al., Ratios of bottom meson branching fractions involving 
J/if mesons and determination of b quark fragmentation fractions, Phys. Rev. D54 
(1996) 6596, arXiv:hep-ex/9607003, 

[32] W. Altmannshofer and D. M. Straub, Implications of b —>■ s measurements, in 50th 
Rencontres de Moriond on EW Interactions and Unified Theories La Thuile, Italy, 
March 14-21, 2015, 2015. arXiv: 1503.06199. 

[33] R. R. Horgan, Z. Liu, S. Meinel, and M. Wingate, Lattice QCD calculation of form 
factors describing the rare decays B —s- K*£ + £~ and B s —> (p£ + £~, Phys. Rev. D89 
(2014) 094501, arXiv:1310.3722, 

[34] R. R. Horgan, Z. Liu, S. Meinel, and M. Wingate, Rare B decays using lattice QCD 
form factors, PoS(LATTICE2014)372 (2015) arXiv: 1501.00367. 

[35] R. R. Horgan, Z. Liu, S. Meinel, and M. Wingate, Calculation of B° —» K*°p + /a~ 
and B —> (fp, + pT observables using form factors from lattice QCD, Phys. Rev. Lett. 
112 (2014) 212003, arXiv:1310.3887. 

[36] A. Ali, E. Lunghi, C. Greub, and G. Hiller, Improved model independent analy¬ 
sis of semileptonic and radiative rare B decays, Phys. Rev. D66 (2002) 034002, 

arXiv:hep-ph/0112300, 

[37] P. Ball and R. Zwicky, B ds —> p,u>,K*,(f decay form factors from light-cone sum 
rules reexamined, Phys. Rev. D71 (2005) 014029, arXiv:hep-ph/0412079, 

[38] S. Descotes-Genon, J. Matias, and J. Virto, Understanding the b —» K*p + p~ anomaly, 
Phys. Rev. D 88 (2013) 074002. 

[39] LHCb collaboration, R. Aaij et al, Precision measurement of CP violation in B ° —> 
J/ifK + K~ decays, Phys. Rev. Lett. 114 (2015) 041801, arXiv: 1411.3104. 

[40] K. De Bruyn et al, Branching ratio measurements of B s decays, Phys. Rev. D86 
(2012) 014027, arXiv:1204.1735, 

[41] G. J. Feldman and R. D. Cousins, A unified approach to the classical statistical 
analysis of small signals, Phys. Rev. D57 (1998) 3873, arXiv:physics/9711021. 

[42] B. Sen, M. Walker, and M. Woodroofe, On the unified method with nuisance parameters, 
Statist. Sinica 19 (2009) 301. 


30 






LHCb collaboration 


R. Aaij 38 , B. Adeva 37 , M. Adinolfi 46 , A. Affolder 52 , Z. Ajaltouni 5 , S. Akar 6 , J. Albrecht 9 , 

F. Alessio 38 , M. Alexander 51 , S. Ali 41 , G. Alkhazov 30 , P. Alvarez Cartelle 53 , A.A. Alves Jr 57 , 

S. Amato 2 , S. Amerio 22 , Y. Arnhis 7 , L. An 3 , L. Anderlini 17,9 , J. Anderson 40 , G. Andreassi 39 , 

M. Andreotti 16 ’-^, J.E. Andrews 58 , R.B. Appleby 54 , O. Aquines Gutierrez 10 , F. Archilli 38 , 

P. d’Argent 11 , A. Artamonov 35 , M. Artuso 59 , E. Aslanides 6 , G. Auriemma 25 ’”, M. Baalouch 5 , 

S. Bachmann 11 , J.J. Back 48 , A. Badalov 36 , C. Baesso 60 , W. Baldini 16,38 , R.J. Barlow 54 , 

C. Barschel 38 , S. Barsuk', W. Barter 38 , V. Batozskaya 28 , V. Battista 39 , A. Bay 39 , L. Beaucourt 4 , 
J. Beddow 51 , F. Bedeschi 23 , I. Bediaga 1 , L.J. Bel 41 , V. Bellee 39 , I. Belyaev 31 , E. Ben-Haim 8 , 

G. Bencivenni 18 , S. Benson 38 , J. Benton 46 , A. Berezhnoy 32 , R. Bernet 40 , A. Bertolin 22 , 

M. -O. Bettler 38 , M. van Beuzekom 41 , A. Bien 11 , S. Bifani 45 , T. Bird 54 , A. Birnkraut 9 , 

A. Bizzeti 17,1 , T. Blake 48 , F. Blanc 39 , J. Blouw 10 , S. Blusk 59 , V. Bocci 25 , A. Bondar 34 , 

N. Bondar 30,38 , W. Bonivento 15 , S. Borghi 54 , M. Borsato 7 , T.J.V. Bowcock 52 , E. Bowen 40 , 

C. Bozzi 16 , S. Braun 11 , D. Brett 54 , M. Britsch 10 , T. Britton 59 , J. Brodzicka 54 , N.H. Brook 46 , 

A. Bursche 40 , J. Buytaert 38 , S. Cadeddu 15 , R. Calabrese 16 ’-^, M. Calvi 20,fc , M. Calvo Gomez 36,p , 
P. Campana 18 , D. Campora Perez 38 , L. Capriotti 54 , A. Carbone 14 ’^, G. Carboni 24,z , 

R. Cardinale 19 ’- 7 , A. Cardini 15 , P. Carniti 20 , L. Carson 50 , K. Carvalho Akiba 2,38 , G. Casse 52 , 

L. Cassina 20,fc , L. Castillo Garcia 38 , M. Cattaneo 38 , Ch. Cauet 9 , G. Cavallero 19 , R. Cenci 23 ’*, 

M. Charles 8 , Ph. Charpentier 38 , M. Chefdeville 4 , S. Chen 54 , S.-F. Cheung 55 , N. Chiapolini 40 , 
M. Chrzaszcz 40 , X. Cid Vidal 38 , G. Ciezarek 41 , P.E.L. Clarke 50 , M. Clemencic 38 , H.V. Cliff 47 , 

J. Closier 38 , V. Coco 38 , J. Cogan 6 , E. Cogneras 5 , V. Cogoni 15,e , L. Cojocariu 29 , G. Collazuol 22 , 
P. Collins 38 , A. Comerma-Montells 11 , A. Contu 15,38 , A. Cook 46 , M. Coombes 46 , S. Coquereau 8 , 
G. Corti 38 , M. Corvo 16 ’^, B. Couturier 38 , G.A. Cowan 50 , D.C. Craik 48 , A. Crocombe 48 , 

M. Cruz Torres 60 , S. Cunliffe 53 , R. Currie 53 , C. D’Ambrosio 38 , E. Dall’Occo 41 , J. Dalseno 46 , 
P.N.Y. David 41 , A. Davis 57 , K. De Bruyn 41 , S. De Capua 54 , M. De Cian 11 , J.M. De Miranda 1 , 

L. De Paula 2 , P. De Simone 18 , C.-T. Dean 51 , D. Decamp 4 , M. Deckenhoff 9 , L. Del Buono 8 , 

N. Deleage 4 , M. Denmier 9 , D. Derkach 55 , O. Deschamps 5 , F. Dettori 38 , B. Dey 21 , A. Di Canto 38 , 
F. Di Ruscio 24 , H. Dijkstra 38 , S. Donleavy 52 , F. Dordei 11 , M. Dorigo 39 , A. Dosil Suarez 37 , 

D. Dossett 48 , A. Dovbnya 43 , K. Dreimanis 52 , L. Dufour 41 , G. Dujany 54 , F. Dupertuis 39 , 

P. Durante 38 , R. Dzhelyadin 35 , A. Dziurda 26 , A. Dzyuba 30 , S. Easo 49,38 , U. Egede 53 , 

V. Egorychev 31 , S. Eidelman 34 , S. Eisenhardt 50 , U. Eitschberger 9 , R. Ekelhof 9 , L. Eklund 51 , 

L El Rifai 5 , Ch. Elsasser 40 , S. Ely 59 , S. Esen 11 , H.M. Evans 47 , T. Evans 55 , A. Falabella 14 , 

C. Farber 38 , C. Farinelli 41 , N. Farley 45 , S. Farry 52 , R. Fay 52 , D. Ferguson 50 , 

V. Fernandez Albor 3 , F. Ferrari 14 , F. Ferreira Rodrigues 1 , M. Ferro-Luzzi 38 , S. Filippov 33 , 

M. Fiore 16,38 ’^, M. Fiorini 16 ’^, M. Firlej 27 , C. Fitzpatrick 39 , T. Fiutowski 27 , K. Fohl 38 , P. Fol 53 , 
M. Fontana 10 , F. Fontanelli 19 ’- 7 , R. Forty 38 , O. Francisco 2 , M. Frank 38 , C. Frei 38 , M. Frosini 17 , 
J. Fu 21 , E. Furfaro 24,1 , A. Gallas Torreira 37 , D. Galli 14,d , S. Gallorini 22,38 , S. Gambetta 50 , 

M. Gandelman 2 , P. Gandini 55 , Y. Gao 3 , J. Garcia Pardihas 37 , J. Garra Tico 47 , L. Garrido 36 , 

D. Gascon 36 , C. Gaspar 38 , R. Gauld 55 , L. Gavardi 9 , G. Gazzoni 5 , A. Geraci 21 ’®, D. Gerick 11 , 

E. Gersabeck 11 , M. Gersabeck 54 , T. Gershon 48 , Ph. Ghez 4 , A. Gianelle 22 , S. Giani 39 , 

V. Gibson 47 , O. G. Girard 39 , L. Giubega 29 , V.V. Gligorov 38 , C. Gobel 60 , D. Golubkov 31 , 

A. Golutvin 53,31 ’ 38 , A. Gomes 1 ’®, C. Gotti 20,fc , M. Grabalosa Gandara 5 , R. Graciani Diaz 36 , 
L.A. Granado Cardoso 38 , E. Grauges 36 , E. Graverini 40 , G. Graziani 17 , A. Grecu 29 , 

E. Greening 55 , S. Gregson 47 , P. Griffith 45 , L. Grillo 11 , O. Griinberg 63 , B. Gui 59 , E. Gushchin 33 , 
Yu. Guz 35 ’ 38 , T. Gys 38 , T. Hadavizadeh 55 , C. Hadjivasiliou 59 , G. Haefeli 39 , C. Haen 38 , 


31 



S.C. Haines 47 , S. Hall 53 , B. Hamilton 58 , X. Han 11 , S. Hansmann-Menzemer 11 , N. Harnew 55 , 
S.T. Harnew 46 , J. Harrison 54 , J. He 38 , T. Head 39 , V. Heijne 41 , K. Hennessy 52 , P. Henrard 5 , 

L. Henry 8 , J.A. Hernando Morata 37 , E. van Herwijnen 38 , M. Hefi 63 , A. Hicheur 2 , D. Hill 55 , 

M. Hoballah 5 , C. Hombach 54 , W. Hulsbergen 41 , T. Humair 53 , N. Hussain 55 , D. Hutchcroft 52 , 

D. Hynds 51 , M. Idzik 27 , P. Ilten 56 , R. Jacobsson 38 , A. Jaeger 11 , J. Jalocha 55 , E. Jans 41 , 

A. Jawahery 58 , F. Jing 3 , M. John 55 , D. Johnson 38 , C.R. Jones 47 , C. Joram 38 , B. Jost 38 , 

N. Jurik 59 , S. Kandybei 43 , W. Kanso 6 , M. Karacson 38 , T.M. Karbach 38 ’’*', S. Karodia 51 , 

M. Kelsey 59 , I.R. Kenyon 45 , M. Kenzie 38 , T. Ketel 42 , B. Khanji 20,38,fc , C. Khurewathanakul 39 , 

S. Klaver 54 , K. Klimaszewski 28 , O. Kochebina 7 , M. Kolpin 11 , I. Komarov 39 , R.F. Koopman, 
P. Koppenburg 41,38 , M. Kozeiha 5 , L. Kravchuk 33 , K. Kreplin 11 , M. Kreps 48 , G. Krocker 11 , 

P. Krokovny 34 , F. Kruse 9 , W. Kucewicz 26,0 , M. Kucharczyk 26 , V. Kudryavtsev 34 , A. 

K. Kuonen 39 , K. Kurek 28 , T. Kvaratskheliya 31 , D. Lacarrere 38 , G. Lafferty 54 , A. Lai 15 , 

D. Lambert 50 , G. Lanfranchi 18 , C. Langenbruch , B. Langhans 38 , T. Latham 48 , C. Lazzeroni 45 , 
R. Le Gac 6 , J. van Leerdam 41 , J.-P. Lees 4 , R. Lefevre 5 , A. Leflat 32,38 , J. Lefrangois 7 , O. Leroy 6 , 

T. Lesiak 26 , B. Leverington 11 , Y. Li', T. Likhomanenko 65,64 , M. Liles 52 , R. Lindner 38 , C. Linn 38 , 

F. Lionetto 40 , B. Liu 15 , X. Liu 3 , D. Loh 48 , S. Lohn 38 , I. Longstaff 51 , J.H. Lopes 2 , 

D. Lucchesi 22,r , M. Lucio Martinez 37 , H. Luo 50 , A. Lupato 22 , E. Luppi 16, Q O. Lupton 55 , 

N. Lusardi 21 , F. Machefert 7 , F. Maciuc 29 , O. Maev 30 , K. Maguire 54 , S. Malde 55 , A. Malinin 64 , 

G. Manca 7 , G. Mancinelli 6 , P. Manning 59 , A. Mapelli 38 , J. Maratas 5 , J.F. Marchand 4 , 

U. Marconi 14 , C. Marin Benito 36 , P. Marino 23,38 ’*, R. Marki 39 , J. Marks 11 , G. Martellotti 25 , 

M. Martin 6 , M. Martinelli 39 , D. Martinez Santos 37 , F. Martinez Vidal 66 , D. Martins Tostes 2 , 

A. Massafferri 1 , R. Matev 38 , A. Mathad 48 , Z. Mathe 38 , C. Matteuzzi 20 , K. Matthieu 11 , 

A. Mauri 40 , B. Maurin 39 , A. Mazurov 45 , M. McCann 53 , J. McCarthy 45 , A. McNab 54 , 

R. McNulty 12 , B. Meadows 57 , F. Meier 9 , M. Meissner 11 , D. Melnychuk 28 , M. Merk 41 , 

D.A. Milanes 62 , M.-N. Minard 4 , D.S. Mitzel 11 , J. Molina Rodriguez 60 , I.A. Monroy 62 , 

S. Monteil 5 , M. Morandin 22 , P. Morawski 27 , A. Morda 6 , M.J. Morello 23,1 , J. Moron 27 , 

A. B. Morris 50 , R. Mountain 59 , F. Muheim 50 , J. Miiller 9 , K. Muller 40 , V. Muller 9 , M. Mussini 14 , 

B. Muster 39 , P. Naik 46 , T. Nakada 39 , R. Nandakumar 49 , A. Nandi 55 , I. Nasteva 2 , 

M. Needham 50 , N. Neri 21 , S. Neubert 11 , N. Neufeld 38 , M. Neuner 11 , A.D. Nguyen 39 , 

T. D. Nguyen 39 , C. Nguyen-Mau 39,9 , V. Niess 5 , R. Niet 9 , N. Nikitin 32 , T. Nikodem 11 , D. Ninci 23 , 
A. Novoselov 35 , D.P. O’Hanlon 48 , A. Oblakowska-Mucha 27 , V. Obraztsov 35 , S. Ogilvy 51 , 

O. Okhrimenko , R. 01deman 15,e , C.J.G. Onderwater 67 , B. Osorio Rodrigues 1 , 

J.M. Otalora Goicochea 2 , A. Otto 38 , P. Owen 53 , A. Oyanguren 66 , A. Palano 13,c , F. Palombo 21, “, 

M. Palutan 18 , J. Panman 38 , A. Papanestis 49 , M. Pappagallo 51 , L.L. Pappalardo 16 ’-', 

C. Pappenheimer 57 , C. Parkes 54 , G. Passaleva 17 , G.D. Patel 52 , M. Patel’ 3 , C. Patrignani 19J , 

A. Pearce 54,49 , A. Pellegrino 41 , G. Penso 25,m , M. Pepe Altarelli 38 , S. Perazzini 14,d , P. Perret 5 , 

L. Pescatore 45 , K. Petridis 46 , A. Petrolini 19 ’- 7 , M. Petruzzo 21 , E. Picatoste Olloqui 36 , 

B. Pietrzyk 4 , T. Pilar 48 , D. Pinci 25 , A. Pistone 19 , A. Piucci 11 , S. Playfer 50 , M. Plo Casasus 37 , 
T. Poikela 38 , F. Polci 8 , A. Poluektov 48,34 , I. Polyakov 31 , E. Polycarpo 2 , A. Popov 35 , 

D. Popov 10,38 , B. Popovici 29 , C. Potterat 2 , E. Price 46 , J.D. Price 52 , J. Prisciandaro 39 , 

A. Pritchard 52 , C. Prouve 46 , V. Pugatch 44 , A. Puig Navarro 39 , G. Punzi 23,s , W. Qian 4 , 

R. Quagliani 7,46 , B. Rachwal 26 , J.H. Rademacker 46 , M. Rama 23 , M.S. Rangel 2 , I. Raniuk 43 , 

N. Rauschmayr 38 , G. Raven 42 , F. Redi 53 , S. Reichert 54 , M.M. Reid 48 , A.C. dos Reis 1 , 

S. Ricciardi 49 , S. Richards 46 , M. Rihl 38 , K. Rinnert 52 , V. Rives Molina 36 , P. Robbe 7,38 , 

A.B. Rodrigues 1 , E. Rodrigues 54 , J.A. Rodriguez Lopez 62 , P. Rodriguez Perez 54 , S. Roiser 38 , 

V. Romanovsky 35 , A. Romero Vidal 37 , J. W. Ronayne 12 , M. Rotondo 22 , J. Rouvinet 39 , T. Ruf 38 , 


32 



H. Ruiz 36 , P. Ruiz Vails 66 , J.J. Saborido Silva 37 , N. Sagidova 30 , P. Sail 51 , B. Saitta 15,e , 

V. Salustino Guimaraes 2 , C. Sanchez Mayordomo 66 , B. Sanmartin Sedes 37 , R. Santacesaria 25 , 

C. Santamarina Rios 37 , M. Santimaria 18 , E. Santovetti 24,i , A. Sarti 18,m , C. Satriano 25,n , 

A. Satta 24 , D.M. Saunders 46 , D. Savrina 31,32 , M. Schiller 38 , H. Schindler 38 , M. Schlupp 9 , 

M. Schmelling 10 , T. Schmelzer 9 , B. Schmidt 38 , O. Schneider 39 , A. Schopper 38 , M. Schubiger 39 , 
M.-H. Scheme 7 , R. Schwemmer 38 , B. Sciascia 18 , A. Sciubba 25,m , A. Semennikov 31 , N. Serra 40 , 

J. Serrano 6 , L. Sestini 22 , P. Seyfert 20 , M. Shapkin 35 , I. Shapoval 16,43 ’^, Y. Shcheglov 30 , 

T. Shears 52 , L. Shekhtman 34 , V. Shevchenko 64 , A. Shires 9 , B.G. Siddi 16 , R. Silva Coutinho 48 , 
G. Sirni 22 , M. Sirendi 47 , N. Skidmore 46 , I. Skillicorn 51 , T. Skwarnicki 59 , E. Smith 55,49 , 

E. Smith 53 , I. T. Smith 50 , J. Smith 47 , M. Smith 54 , H. Snoek 41 , M.D. Sokoloff 57,38 , F.J.P. Soler 51 , 

F. Soornro 39 , D. Souza 46 , B. Souza De Paula 2 , B. Spaan 9 , P. Spradlin 51 , S. Sridharan 38 , 

F. Stagni 38 , M. Stahl 1 11 , S. Stahl 38 , O. Steinkamp 40 , O. Stenyakin 35 , F. Sterpka 59 , 

S. Stevenson 55 , S. Stoica 29 , S. Stone 59 , B. Storaci 40 , S. Stracka 23,4 , M. Straticiuc 29 , 

U. Straumann 40 , L. Sun 57 , W. Sutcliffe 53 , K. Swientek 27 , S. Swientek 9 , V. Syropoulos 42 , 

M. Szczekowski 28 , P. Szczypka 39,38 , T. Szumlak 27 , S. T’Jampens 4 , A. Tayduganov 6 , 

T. Tekampe 9 , M. Teklishyn 7 , G. Tellarini 16 ^, F. Teubert 38 , C. Thomas 55 , E. Thomas 38 , 

J. van Tilburg 41 , V. Tisserand 4 , M. Tobin 39 , J. Todd 57 , S. Tolk 42 , L. Tomassetti 16 ’^, 

D. Tonelli 38 , S. Topp-Joergensen 55 , N. Torr 55 , E. Tournefier 4 , S. Tourneur 39 , K. Trabelsi 39 , 
M.T. Tran 39 , M. Tresch 40 , A. Trisovic 38 , A. Tsaregorodtsev 6 , P. Tsopelas 41 , N. Tuning 41,38 , 

A. Ukleja 28 , A. Ustyuzhanin 65,64 , U. Uwer 11 , C. Vacca 15,e , V. Vagnoni 14 , G. Valenti 14 , 

A. Vallier', R. Vazquez Gomez 18 , P. Vazquez Regueiro 37 , C. Vazquez Sierra 37 , S. Vecchi 16 , 

J.J. Velthuis 46 , M. Veltri 17,/l , G. Veneziano 39 , M. Vesterinen 11 , B. Viaud 7 , D. Vieira 2 , 

M. Vieites Diaz 37 , X. Vilasis-Car dona 36 ,p , A. Vollhardt 40 , D. Volyanskyy 10 , D. Voong 46 , 

A. Vorobyev 30 , V. Vorobyev 34 , C. Vofi 63 , J.A. de Vries 41 , R. Waldi 63 , C. Wallace 48 , R. Wallace 12 , 
J. Walsh 23 , S. Wandernoth 11 , J. Wang 59 , D.R. Ward 47 , N.K. Watson 45 , D. Websdale 53 , 

A. Weiden 40 , M. Whitehead 48 , G. Wilkinson 55,38 , M. Wilkinson 59 , M. Williams 38 , 

M.P. Williams 45 , M. Williams 56 , T. Williams 45 , F.F. Wilson 49 , J. Wimberley 58 , J. Wishahi 9 , 

W. Wislicki 28 , M. Witek 26 , G. Wormser 7 , S.A. Wotton 47 , S. Wright 47 , K. Wyllie 38 , Y. Xie 61 , 

Z. Xu 39 , Z. Yang 3 , J. Yu 61 , X. Yuan 34 , O. Yushchenko 35 , M. Zangoli 14 , M. Zavertyaev 10,5 , 

L. Zhang 3 , Y. Zhang 3 , A. Zhelezov 11 , A. Zhokhov 31 , L. Zhong 3 , S. Zucchelli 14 . 

1 Centro Brasileiro de Pesquisas Fisicas (CBPF), Rio de Janeiro, Brazil 

2 Universidade Federal do Rio de Janeiro (UFRJ), Rio de Janeiro, Brazil 

3 Center for High Energy Physics, Tsinghua University, Beijing, China 

4 LAPP, Universite Savoie Mont-Blanc, CNRS/IN2P3, Annecy-Le-Vieux, France 
5 Clermont Universite, Universite Blaise Pascal, CNRS/IN2P3, LPC, Clermont-Ferrand, France 
6 CPPM, Aix-Marseille Universite, CNRS/IN2P3, Marseille, France 
7 LAL, Universite Paris-Sud, CNRS/IN2P3, Orsay, France 

8 LPNHE, Universite Pierre et Marie Curie, Universite Paris Diderot, CNRS/IN2P3, Paris, France 

9 Fakultat Physik, Technische Universitat Dortmund, Dortmund, Germany 

10 Max-Planck-Institut fiir Kernphysik (MPIK), Heidelberg, Germany 

11 Physikalisches Institut, Ruprecht-Karls-Universitat Heidelberg, Heidelberg, Germany 

12 School of Physics, University College Dublin, Dublin, Ireland 

13 Sezione INFN di Bari, Bari, Italy 

1A Sezione INFN di Bologna, Bologna, Italy 

15 Sezione INFN di Cagliari, Cagliari, Italy 

16 Sezione INFN di Ferrara, Ferrara, Italy 

17 Sezione INFN di Firenze, Firenze, Italy 


33 



18 Laboratori Nazionali dell’INFN di Frascati, Frascati, Italy 

19 Sezione INFN di Genova, Genova, Italy 

20 Sezione INFN di Milano Bicocca, Milano, Italy 

21 Sezione INFN di Milano, Milano, Italy 

22 Sezione INFN di Padova, Padova, Italy 

23 Sezione INFN di Pisa, Pisa, Italy 

24 Sezione INFN di Roma Tor Vergata, Roma, Italy 

25 Sezione INFN di Roma La Sapienza, Roma, Italy 

26 Henryk Niewodniczanski Institute of Nuclear Physics Polish Academy of Sciences, Krakow, Poland 
2 ‘ AGH - University of Science and Technology, Faculty of Physics and Applied Computer Science, 
Krakow, Poland 

28 National Center for Nuclear Research (NCBJ), Warsaw, Poland 

29 Horia Hulubei National Institute of Physics and Nuclear Engineering, Bucharest-Magurele, Romania 

30 Petersburg Nuclear Physics Institute (PNPI), Gatchina, Russia 

31 Institute of Theoretical and Experimental Physics (ITEP), Moscow, Russia 

32 Institute of Nuclear Physics, Moscow State University (SINP MSU), Moscow, Russia 
33 Institute for Nuclear Research of the Russian Academy of Sciences (INR RAN), Moscow, Russia 
34 Budker Institute of Nuclear Physics (SB RAS) and Novosibirsk State University, Novosibirsk, Russia 

35 Institute for High Energy Physics (IHEP), Protvino, Russia 

36 Universitat de Barcelona, Barcelona, Spain 

37 Universidad de Santiago de Compostela, Santiago de Compostela, Spain 

38 European Organization for Nuclear Research (CERN), Geneva, Switzerland 
39 Ecole Polytechnique Federate de Lausanne (EPFL), Lausanne, Switzerland 

40 Physik-Institut, Universitat Zurich, Zurich, Switzerland 

41 Nikhef National Institute for Subatomic Physics, Amsterdam, The Netherlands 

42 Nikhef National Institute for Subatomic Physics and VU University Amsterdam, Amsterdam, The 
Netherlands 

43 NSC Kharkiv Institute of Physics and Technology (NSC KIPT), Kharkiv, Ukraine 

44 Institute for Nuclear Research of the National Academy of Sciences (KINR), Kyiv, Ukraine 

45 University of Birmingham, Birmingham, United Kingdom 

46 H.H. Wills Physics Laboratory, University of Bristol, Bristol, United Kingdom 

47 Cavendish Laboratory, University of Cambridge, Cambridge, United Kingdom 

48 Department of Physics, University of Warwick, Coventry, United Kingdom 

49 STFC Rutherford Appleton Laboratory, Didcot, United Kingdom 

50 School of Physics and Astronomy, University of Edinburgh, Edinburgh, United Kingdom 
51 School of Physics and Astronomy, University of Glasgow, Glasgow, United Kingdom 
° 2 Oliver Lodge Laboratory, University of Liverpool, Liverpool, United Kingdom 
53 Imperial College London, London, United Kingdom 

54 School of Physics and Astronomy, University of Manchester, Manchester, United Kingdom 
55 Department of Physics, University of Oxford, Oxford, United Kingdom 
56 Massachusetts Institute of Technology, Cambridge, MA, United States 

57 University of Cincinnati, Cincinnati, OH, United States 

58 University of Maryland, College Park, MD, United States 

59 Syracuse University, Syracuse, NY, United States 

60 Pontificia Universidade Catolica do Rio de Janeiro (PUC-Rio), Rio de Janeiro, Brazil, associated to 2 

61 Institute of Particle Physics, Central China Normal University, Wuhan, Hubei, China, associated to 3 

62 Departamento de Fisica , Universidad Nacional de Colombia, Bogota, Colombia, associated to 8 
63 Institut fur Physik, Universitat Rostock, Rostock, Germany, associated to 11 

"National Research Centre Kurchatov Institute, Moscow, Russia, associated to 31 
65 Yandex School of Data Analysis, Moscow, Russia, associated to 31 

"institute de Fisica Corpuscular (IFIC), Universitat de Valencia-CSIC, Valencia, Spain, associated to 36 
67 Van Swinderen Institute, University of Groningen, Groningen, The Netherlands, associated to 41 


34 



a Universidade Federal do Triangulo Mineiro (UFTM), Uberaba-MG, Brazil 

b P.N. Lebedev Physical Institute, Russian Academy of Science (LPI RAS), Moscow, Russia 

c Universita di Bari, Bari, Italy 

d Universita di Bologna, Bologna, Italy 

e Universita di Cagliari, Cagliari, Italy 

f Universita di Ferrara, Ferrara, Italy 

9 Universita di Firenze, Firenze, Italy 

h Universita di Urbino, Urbino, Italy 

1 Universita di Modena e Reggio Emilia, Modena, Italy 

0 Universita di Genova, Genova, Italy 

k Universita di Milano Bicocca, Milano, Italy 

1 Universita di Roma Tor Vergata, Roma, Italy 

m Universita di Roma La Sapienza, Roma, Italy 

n Universita della Basilicata, Potenza, Italy 

°AGH - University of Science and Technology, Faculty of Computer Science, Electronics and 
Telecommunications, Krakow, Poland 

p LIFAELS, La Salle, Universitat Ramon Hull, Barcelona, Spain 
q Hanoi University of Science, Hanoi, Viet Nam 
r Universita di Padova, Padova, Italy 
s Universita di Pisa, Pisa, Italy 
4 Scuola Normale Superiore, Pisa, Italy 
“ Universita degli Studi di Milano, Milano, Italy 
v Politecnico di Milano, Milano, Italy 

t Deceased 


35 



